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ABSTRACT 

It has been known for over 30 years that Galactic globular clusters (GCs) are overabundant by 
orders of magnitude in bright X-ray sources per unit mass relative to the disk population. Recently 
a quantitative understanding of this phenomenon has developed, with a clear correlation between 
the number of X-ray sources in a cluster. Nx, and the cluster's encounter frequency, F, becoming 
apparent. We derive a refined version of F that incorporates the finite lifetime of X-ray sources and 
the dynamical evolution of clusters. With it we find we are able to explain the few clusters that lie 
off the iVx-F correlation, and resolve the discrepancy between observed GC core radii and the values 
predicted by theory. Our results suggest that most GCs are still in the process of core contraction 
and have not yet reached the thermal equilibrium phase driven by binary scattering interactions. 
Subject headings: globular clusters: general — methods: numerical — stellar dynamics 



1. X-RAY SOURCES AND CLUSTER DYNAMICS 



It was realized more than 30 years ago that Galactic 
globular clusters (GCs) are overabundant by orders of 
magnitude in bright X-r ay sources per unit mas s relative 
to the disk population ()Clarklll97a iKatd 1X9751 ). It was 
quickly understood that strong dynamical scattering in- 
teractions of binaries in the dense clus ter cores should 
be re sponsible for this overabundance (jVerbunt fc HutI 
119871) . With the advances in X-ray astronomy made pos- 
sible by observatories such as Chandra, the relationship 
between X-ray sources and core clus ter dynamics has re - 
cently been quantitatively studied. iPoolev et al.l (120031 ) 
performed Chandra observations of many Galactic GCs 
down to a limiting luminosity of 4 x lO'^" erg/s in the 
0.5-6 keV range (which includes low-mass X-ray bina- 
ries [LMXBs] in outburst and quiescence, cataclysmic 
variables [CVs] , millisecond pulsars [MSPs] , and magnet- 
ically active main sequence binaries [ABs]), and looked 
for correlations between the number of X-ray sources in 
each cluster and properties of the cluster itself. They 
found the strongest correlation with the "encounter fre- 
quency" F, a rough estimate of the current dynamical en- 
counte r rat e in the cluster . Mor e recentlv. iHeinke et al.l 
(|2003f l and IPoolev fc Rv^ (|2006[ ) have isolated the qui- 
escent LMXBs (qLMXBs) and CVs, respectively, from 
the X-ray source populations, and have shown that their 
numbers are indeed consistent with dynamical formation. 

These results represent quantitative, empirical evi- 
dence that dynamical encounters are responsible for the 
formation of X-ray sources in clusters. However, they 
suffer from at least a few drawbacks. First, the corre- 
lation between the number of X-ray sources, Nx, and 
the encounter frequency appears to be sub-linear, with 
Nx oc F°-^ ^^°-3^, although fo r LMXBs the exponent is 
0.97 ±0.5 (jPoolev et al.ll2003[ ). Second, there are three 
clusters for which Nx is significan tly larger than pre - 
dicted by F. In the original work of IPoolev et al.l (j2003l ) 
it was already clear that NGC 6397 has an Nx that is 
~ 5 times larger than predicted by the Nx-T correla- 
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Fig. 1. — Number of observed cluster X-ray sources with Lx 
4 X 10^'^ erg/s for several Galactic GCs versus the e ncount e r rate 
r. The power-law fit and data poin ts are from Poolcy JTaD li2003) 
with the exception of NGC 7099 IILueger et al.n2007l) and Ter 1 
l|Cackett et al.ll2006l 'l. The Nx error bars for Ter 1, NGC 6397, 
and NGC 7099 represent source counting noise and background 
source uncertainty, but for the remaining clusters represent only 
background uncertainty. 



tion. Recent observations show that Nx is factor of ^ 2 
times that predicted by F for NGC 7099 (iLugger et al 
2007D . and a factor of ~ 20 for Ter 1 (jCackett et al 



20m . The common thread among these three clusters 
is that they are o bservationally "core -collapsed," while 
all oth ers in the IPoolev et al.l (|2003f ) sample are not 
(jGrindla v 200^). (A possible exception i s NGC 6752 
whose collapsed core st atus is debated (jFerraro et al.l 
120031 : iLugger et a"nil995[ ).') A cluster is observationally 
termed core-collapsed if its surface brightness profile is 
consistent with a cusp at the limit of resolution, making 
it more difficult to measure the core radius. As described 
below, the collapsed core status of a cluster is linked to 
its dynamical state, implying that cluster evolution com- 
plicates the Nx-T correlation. 

2. UNDERSTANDING CLUSTER CORE RADII 

The evolution of a GC, being a bound self-gravitating 
system, is very similar to the evolution of a star, and 
comprises three main phases. In the "core-contraction 
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phase," the first phase of evolution, tlie cluster's core 
contracts on a relaxation timescale, much like a pre- 
main sequence star. Once the core density becomes large 
enough for binary stars to begin strongly interacting dy- 
namically, and thus generating energy via super-elastic 
encounters, the cluster settles into the "binary-burning 
phase," analogous to the main sequence in stars. In this 
phase the dynamical pro perties of the cluster core remain 
roughly constant (e.g., iGao et al.l llQQlt iHeggie fc HutI 
120031: iFregeau fc Rasiol[2007[ ). Once the binary popula- 
tion is exhausted in the cluster core, it will collapse'^ via 
the gravothermal instability, leading to extremely high 
central densities. Deep in collapse, an energy producing 
event, such as an interaction of a dynamically- formed 
binary, will reverse the collapse, causing the core to re- 
bound and enter the "gravothermal oscillation phase," 
in which the core c ontinues to collapse and rebound 
(jHeggie fc HutI [20031) . For a graphical representation of 
th e three ma i n pha ses of cluste r evolution, se e Figu re 1 
of lGao et al.l (fl99ll). Figure 5 of'Frcg eau etall (|2003[ ). or 
Figure 29.1 of iHeggie fc mrt (2003). 

Since there is strong observational evidence that Galac- 
tic GCs were b orn with significant binary fractions 
(|Hut et al.lll992f ). and since the binary-burning phase is 
the longest-lived p hase of cluster eyolution (p erhaps tens 
of Hubble times: iGao et all ll99lUHeggie fc H ut 200l; 
IFregeau fc Rasiol(2007l : rTrenti et al.ll2007r ), it is widelv be- 
lieved that most clusters observed today should be in this 
phase. Early approximate calculations suggested that 
the ratio of core to half-mass radius, rc/r/j, in the binary 
burning phase, is broadly consistent with observations of 
the ~ 80% of Galactic GCs that are not observatio nally 
core-collapsed (|Gao et al.lll991l : IFregeau et all2003f ). Re- 
cently, however, more accurate simulations have shown 
that the early calc ulations overe stimate rdrh by a fac- 
tor o f 10 or more (jHeggie et al.l 2006; F regeau fc Rasiol 
120071 ). These latest results are quite difficult to ignore, 
since they represent the concordance of two completely 
independent cluster evolution codes — one direct A^-body, 
with minimal approximations and a natural inclusion of 
binary interactions; the other the approximate Hcnon 
Monte Carlo method, with direct few-body integration 
of binary interactions. The values of rdrh in the binary- 
burning phase agree quite well between the two codes, 
and fur thermore agree quite well with semi-analytical 
theory (jVesperini fc Cherno3ll994f ). The result is that 
now only the core-collapsed clusters agree with the pre- 
dicted values of r^/rh in the binary-burning phase, im- 
plying that if most clusters arc in this phase some other 
energy generation mechanism is responsible for the mea- 
sured core sizes. Several suggestions have been put forth 
for the energy source, including, most notably, cent ral 
intermediate-mass black holes (IMBHs) (|Trentill2006[) . 

One key feature of the evolution of GCs that has es- 
caped careful attention, though, is the timescale of the 
initial phase of core contraction. Many numerical simula- 
tions have shown that it is of order ^10 initial half-mass 

^ Note that due to clashing naming conventions between observa- 
tion and theory, this core collapse docs not necessarily correspond 
to the observational definition, as discussed in the text. We there- 
fore suggest that some non-dynamical term be used in place of the 
observational "core-collapsed" designation, or at least that the dis- 
tinction between the observational and theoretical uses of the term 
be noted when it is used. 



relaxation times ( Gao et al.l [l99ll iGiersz fc Spurzeml 



120031 IHeggie et al.l 120061 : IFregeau fc Rasiol 1200 7( )! Since 
the current half-mass relaxation time for most clusters is 
~ 1 Gyr and was likely much longer in the past, the core 
contractio n phase may easily last longer than a Hubble 
time (e.g.. lHurlevll2007l ) . In fact, recent iV-body simula- 
tions have shown that for a range of initial binary frac- 
tions, the evolution of Vc/rh over a Hubble time is a sim- 
ple core contraction w ith no evidenc e of a binary-burning 
phase being reached (jHurlevI l2007f ). Thus the solution 
to the discrepancy between theory and observations in 
Tc/rh could be the most mundane one, namely that core- 
collapsed clusters are in the binary-burning phase while 
the rest are still undergoing core contraction. 

3. A REFINED V 

First introduced by IVerbunt fc HutI ()1987[ ). the en- 
counter frequency F is an estimate of the current dy- 
namical interaction rate in the cluster, which is assumed 
to be proportional to the current number of observable 
X-ray sources. We refine this predictive quantity by in- 
cluding, among other factors, the recent history of the 
evolution of the core properties. We start with the gen- 
eral form of the inter action rate, then specialize to the 
IVerbunt fc H^ (IT987I ) F, and our new version. 

The interaction rate between two species of objects can 
be written generally as 



dNin 

dt 



nin2CTi2 |vi2|/(vi2)(i^vi2(i^r, (1) 



where is the number density of species i, CTi2 is the 
interaction cross section between the two species, Vi2 is 
their relative speed with / its distribution function, and 
the integral is carried out over relative velocity space and 
volume. For the dynamical creation of LMXBs, species 1 
represents stellar binaries, while species 2 represents neu- 
tron stars. The dynamical formation of CVs and other 
low-luminosity X-ray sources is rather more complicated, 
so species 1 and 2 generally repres ent single and binar y 
star systems (see, e.g.. Figure 7 of llvanova et al.ll2006l) . 
Typically, the integral in eq. ([T]) is approximated as 



dt 



« Plrl/va 



(2) 



where is the core mass density, v„ is the core velocity 
dispersion, the integral has been approximated by the 
core value, the gravitational-focusing dominated cross 
section has been used, and Pc,i/ Pc,2 is assumed to be con- 
stant for all clusters. Additionally, it should be pointed 
out that only the proportionality in eq. ([T|) has been pre- 
served since several factors (some of which are not con- 
stant among clusters) have been dropped. More accurate 
approximations of eq. ([T]) have been used, including nu - 
merical integrals over cluster models (jPoolev et al.l20d3l ). 
but all arc estimates of the current interaction rate. 

Dynamically formed X-ray binaries (XRBs) are known 
to have finite detectable lifetimes. For LMXBs, this life- 
time varies from ~ 10^-10^ yr for red giant donors, to 
~ 1 Gyr for main- sequence companion s, to a few Gyr 
for ultracompacts (jlvanova et al.l [20071 ). For CVs, this 
lifetime is ~ 1 Gyr (N. Ivanova, priv. comm.). An ad- 
ditional complication is that dynamically formed XRBs 
do not necessarily turn on as X-ray sources immediately 
following a strong interaction. In fact, the interaction 
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that places a binary on the path to becoming an observ- 
able X-ray source ty pically occurs several Gy r before it 
becomes detectable (jlvanova et al.|[2006ll200lD . 

Since the XRBs we see now formed a few Gyr or more 
ago, and since the recent dynamical history of clusters 
may have been quite variable, it is clear that the current 
number of observable sources should be proportional to 
the interaction rate integrated over time. We thus write 



phase to the same cluster in core contraction is then 



nin2CTi2 |vi2| /(vi2)d^vi2 (frdt. (3) 



We perform the time integration over an interval t^, the 
typical detectable lifetime of an XRB, to a time tg, the 
typical timescale between strong interaction and observ- 
ability, in the past. We leave tx and ti as parameters, 
but take ~ 1 Gyr and = 3 Gyr as canonical values. 
We further simplify the integral by writing 



fbfconlastiongV2vcr—rldt, (4) 

to-tf-tx 



with fb the core binary fraction, fco the core compact 
object fraction, ric the core number density, Vc the core 
radius, Va the 1-D core velocity dispersion, and iTstrong 
the cross section for a strong interaction between a binary 
and a single star. The factor of -\/2 is from taking the 
difference of two Maxwellian velocity distributions. The 
cross section is (Tstrong ~ Tia2GM / {\/2vaY , where a is the 
typical semi-major axis of a binary, M — 3m is the total 
mass of the binary-single system, and m is the typical 
stellar mass. Taking /& and fco to be constant over the 
time integral and substitut ing the definition of the core 
radius, rl = 9vl/ATrGmnc (|Heggie fc Hutll2003f) . yields 



-^int — fbfco 



81%/2a 
4Gm 



dt. 



(5) 



tQ—ti—tx 



This expression is similar to eq. but with the core 
properties integrated over time. We have kept all nu- 
merical factors for the sake of completeness — only the 
proportionality represented by this equation is needed for 
what follows. By substituting in for the time evolution of 
the core quantities in the integrand, one can use this ex- 
pression to differentiate among the three different phases 
of cluster evolution. We exclude the gravothermal oscil- 
lation phase, since in this phase a clust er should have 
no m ore than a few binaries in its core (jHeggie fc HutI 
l2003l ). and all clusters considered here with measured 
core binary fractions show evidence for many more bina- 
ries than this (note that a recent estima te puts the core 
binar y fraction of NGC 6397 at 15± l%: [Da^is fc Riched 
[2001 . 

In the binary-burning phase the core radius and ve- 
locity dispersion are constant, so the integral is easy 
to evaluate. For the core contraction phase, we adopt 
the time evolut ion of the cor e radius shown in the A^- 
body models of iHurlevI (|2007t ). approximated as rdt) ~ 
?'c,o(10 — 9t/to), with rc,o the current core radius and 
to the current cluster age. Since there is essentially no 
core energy support in the core contraction phase, for 
the relationship between core radius and central velocity 
dispersion w e adopt the self-similar co llapse model, with 
vl oc r-°-2i (jBinnev fc Tremainel[T987l ) . The ratio of the 
number of interactions for a cluster in the binary-burning 
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2.835 



(to 



to 



0.315 



(toH 



0.315 ' 



(6) 



.,+9tf J \to+9te+9t^ 

which has a minimum of 2.0 and a maximum of 17.8 
in the range = 10~^-3Gyr, tg ~ 1-10 Gyr, for 
tQ = 13 Gyr. For the canonical values of = 1 Gyr 
and = 3 Gyr with to — 13 Gyr, the value is 5.0. Since 
the number of X-ray sources should scale roughly linearly 
with the number of interactions, this suggests that if a 
cluster is in the binary-burning phase (and has been for 
a time tg + tx to the present), it should have ~ 5 times as 
many X-ray sources than it would if it were in the core 
contraction phase. (If the iPoolev et al.l (j2003i) exponent 
of 0.74 is adopted this factor is 3.3.) Interestingly, the 
three clusters in Figure [1] that are observed to be core 
collapsed are the three that have a significantly larger 
Nx than predicted by the standard Nx-^ correlation, 
by a factor of 2 to ~ 20. NGC 6397 has an overabun- 
dance factor of 2, NGC 7099 a factor of ~ 5, and Ter 
1 a factor of ~ 20. However, we note that the estimate 
of F for Ter 1 is more uncertain than the rest since it 
comes directly from eq. ([2]) and not an integral over the 
cluster profile, and since the velocity dispersion is only 
estimated as it has not been measured observationally. 
The similarity between the observed overabundances and 
those predicted by our simple revision of F suggests that 
the observationally core-collapsed clusters are indeed in 
the binary-burning phase, while the rest are still in the 
process of core contraction. 

A natural objection to this conclusion is that clus- 
ter metallicity may explain the X-ray source overabun- 
dance. Studies of bright LMXBs associated with star 
clusters in external galaxies have shown a strong corre- 
lation between cluster metallicity and LMXB incidence, 
with metal-rich clusters being much more likely to har- 
bor an LMXB (by a fact or of 3 or more; iJordan et al.l 
l2004HSivakoff et al.l[200l . However, this result appears 
to hold only for bright LMXBs. When looking at the 
low-luminosity LMXB and CV populations in our Galaxy 
(which comprise the majority of sources in Figure [T]), a 
correlation between source i ncidence and cluste r metal- 
licity is not clearly apparent (jHeinke et al.ll2006l ) . In any 
case, the mctallicities of the three overabundant clusters 
are not significantly larger than those of the other clus- 
ters in FigurelH with [Fe/H] values of -1.30, -1.95, -2.12 
for Ter 1, NGC 6397, and NGC 7099, respectively, with 
the rest of the clusters ranging in value from -0.34 (NGC 
6440) to -1.75 (NGC 6093) (Harris 1996). 

Another possibility is that the XRBs in the low-F clus- 
ters may be mainly primordial, in which case their num- 
ber should scale with cluster mass. However, were this 
the case, a look at cluster absolute magnitudes shows 
that NGC 6121 would have roughly as many sources as 
NGC 7099 and more than NGC 6397 and Ter 1. an d NGC 
6366 would have more than Ter 1 (jHarrisI [19961 ). Note 
that the 2-5 sources detected in the low-F cluster NGC 
288 (not included i n this study), with F « 0.53 in the 
iPoolev et al. (l2003f) normalization, are likely primordial 
(|Kong et al.ll2006[ ). 



4. DISCUSSION 

This letter presents the confluence of three suggestive 
observational and theoretical results into a self-consistent 
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picture. The first is that of tlic clusters that have been 
observed sufficiently to determine their XRB population, 
the three that are core-collapsed are the same three that 
have a significant X-ray source overabundance (of a fac- 
tor of ~ 2 to ^ 20). The second is the semi-analytical 
result derived in this letter that a cluster in the binary 
burning phase for the last few Gyr should have ~ 5 times 
more dynamically formed X-ray sources than if it were in 
the core contraction phase for the same time. The third is 
the recently confirmed discrepancy between observations 
and theory for the core radii of Galactic GCs, which sug- 
gests that only the observationally core-collapsed clusters 
are in the binary-burning phase. In light of these facts, 
the conclusion that seems strongly suggested is that most 
Galactic GCs are currently still in the core contraction 
phase, while only the 20% of clusters that are core- 
collapsed are in the binary burning phase. This goes 
counter to the widely held belief that most clusters are 
currently in the binary burning phase, and complicates 
the many existing studies that have assumed cluster core 
properties that are constant with time. 

The implications of this result are manifold. There 
are many studies of the dynamical production of in- 
teresting source populations in clusters which assume 
core properties that are constant with time. These in- 
clude predictions of the formation of blue stragglers (e.g., 
iMapelh et al.l I2004D. the ev olution of the core binary 
fracti on lllvanova et al.ll2005D. and t idal-capture binaries 
(e.g.. lDi Stefano fc Rappaportlll994l ). among others. Re- 
vising the results may be as simple as scaling predicted 
source numbers, but may not be so simple for other quan- 
tities. 

Studies of GC evolution have shown that clusters start- 
ing from very different initial conditions evolve toward 
a common range of values in many observable struc- 
tural parameters in the binary-burning phase, including 
the concentration and ratio of core to half- mass radius 
(iFregeau et al.ll2003tlHeggie et al.l2006HFregeau fc Rasiol 
|2007() . Since most clusters may not be in the binary- 
burning phase after all, their observed properties are 
likely to be more strongly correlated with their initial 
conditions. This makes modeling of clusters a bit more 
complicated, but on the other hand allows one to more 



readily deduce something about the initial properties of 
clusters. 

Perhaps anticlimactically, our results suggest that the 
alternative energy sources recently proposed for sup- 
porting GC cores are not required. These inclu de the 
sugge stion of IMBHs in many Galactic GCs (jTrentil 
|2006[) . enhanced stellar mass lo ss from stellar evolutio n 
of physical collision products (|Chatteriee et al.l l2007l ). 
mass segregation of compact remnants in young clusters 
(jMerritt et al.|[200^ . or evaporation of t he stellar-mass 
black hole subsystem in young clusters (jMackev et al.l 
[2001 . 

Although the picture painted in this letter is a sugges- 
tive one, there are still several caveats and limitations 
to our analysis. In the derivation of our refined F we 
assume that the core binary fraction and abundance of 
compact objects are constant over the time of integra- 
tion. Neither is strictly true, although we expect they 
will not vary enough to significantly change the over- 
abundance value we derive. Wc have used only one pos- 
sible expression for the evolution of the core radius in 
core contraction. While we expect the general behavior 
to be very similar to what we have assumed here, more 
work needs to be done to determine if it is universal. Ad- 
ditionally, our analysis ignores the effect of Galactic tidal 
stripping on cluster mass, which would make a cluster ap- 
pear overabun dant in X-ray sourc es, and may be relevant 
for NGC 6397 (jPoolev et al.ll2003[) . On the observational 
side, there are some uncertainties in evaluating F, which 
is dependent on quantities that are somewhat difficult to 
measure for core-collapsed clusters. 
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